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Abstract 



The number state method is used to study sohton bands for three an- 
harmonic quantum lattices: i) The discrete nonhnear Schrodinger equation, 
ii) The Ablowitz-Ladik system, and iii) A fermionic polaron modeL Each of 
these systems is assumed to have /-fold translational symmetry in one spatial 
dimension, where / is the number of freedoms (lattice points). At the second 
quantum level (n = 2) we calculate exact eigenfunctions and energies of pure 
quantum states, from which we determine binding energy (-Eb), effective mass 
(m*) and maximum group velocity CKn) of the soliton bands as functions of 
the anharmonicity in the limit / — » cx). For arbitrary values of n we have 
asymptotic expressions for E^,, m* , and Vm as functions of the anharmonic- 
ity in the limits of large and small anharmonicity. Using these expressions 
we discuss and describe wave packets of pure eigenstates that correspond to 
classical solitons. 



1 Introduction 

Throughout the development of modern nonlinear dynamics the investigation of 
lattices has played a significant role. Often the motivation for such studies is that 
molecular crystals are lattices, and in these applications quantum effects cannot be 
ignored. Typical experiments — such as infra-red absorption, Raman scattering, and 
neutron diffraction — deal with line spectra, where every line corresponds to a pair 
of quantum states, each with a particular number of quanta. At larger quantum 
numbers — approaching the correspondence limit — one is often interested in knowing 
how quantum corrections alter the results of classical calculations. Our aim in this 
paper is to present some results of an exact theory of lattice solitons in the quantum 
regime. To this end we consider the following specific models, which are defined 
on one- dimensional lattices of / freedoms (or lattice sites) with periodic boundary 
conditions and with 7 being the ratio of anharmonicity to nearest neighbor hopping 
energy. 



(i) The quantum discrete nonlinear Schrodinger (QDNLS) equation 



This system arises in the study of molecular vibrations in one-dimensional chains 
such as benzene and certain molecular crystals ||^, ^ ^ . With site energies 

scaled out through a gauge transformation the reduced Hamiltonian operator is 

/ 

Hi = -J2 [%+^ + + ' (1-1) 

where bj and bj are standard bosonic raising and lowering operators satisfying the 
commutation relations [bi,bj] — — 0, = <5,;j at each freedom. 



(ii) The quantum Ablowitz-Ladik (QAL) equation 



The Ablowitz-Ladik equation is of interest because the corresponding classical 
system is integrable via the inverse scattering method Q and the quantum system 
is a simple example of a g-boson model p^ . The reduced Hamiltonian operator is 

/ 



'1+1 + bj-i) 



(1.2) 



1 



where and bj are operators satisfying "g-deformed" commutation relations 



[bl,bi] = [b„b,] = 0, [b„bl] = (l + lb]b,) 5, 



(iii) A fermionic polaron (FP) model 

This model describes the dynamics of electrons in a one-dimensional crystal 
and is related to the XXZ spin chain model . The reduced Hamiltonian operator 
is 



/ 

i=i 



+ + 7a]aja]+iaj+iJ , (1-3) 



where Oj and aj are standard fermion raising and lowering operators satisfying the 
anticommutation relations {a^, aj} = {a|, aj} = 0, {ci, aj} = Sij. 

These systems are chosen to exhibit three different types of fundamental quanta: 
i) standard bosons, ii) g-deformed bosons, and iii) standard fermions. A classical 
analog (or correspondence limit) exists only in the first two examples. 

Since each of these three models is assumed to satisfy periodic boundary condi- 
tions, the Hamiltonians are invariant under the action of the translation operator 
with eigenvalue exp (ifc) , where k is the crystal momentum. For a fixed number 
of quanta (n) the energy eigenstates for a particular value of k are separated into 
bands, and the band of lowest energy is of the form [| 

E = En{k), (1.4) 

which is of particular interest for two reasons: 

• Within this band there is but a single eigenstate for each value of k. 

• The n quanta of the eigenfunctions for the band are located more closely 
together than in the bands of higher energy. 



We call this lowest band the soliton band and use Equation (1.4) to compute the 
binding energy, effective mass, and maximum group velocity of the corresponding 
quantum soliton. 

In the following section we sketch our procedure for computing En{k) as a func- 
tion of the ratio of anharmonicity to nearest neighbor hopping energy (7) , and in 
Section 3 we give complete descriptions of the bands of our three models for the 
second quantum level (n = 2). Asymptotic expressions for i?„(fc) at small and large 
values of 7 are derived in Section 4. In Section 5 we use these results to construct 
wave packets — over n and k — of pure eigenstates that correspond to classical soli- 
tons, and some conclusions are presented in Section 6. Throughout the discussion 
we assume units of energy and time for which h — 1. 

We have chosen to call our method of analysis the number state method, to 
distinguish it from the Quantum Inverse Scattering Method ||l^, |l^ . The advantages 
and disadvantages of each of these two methods is discussed at some length in . 
In all cases where both techniques have been used they give (as expected) the same 
results. However it is important to note that the first model we consider (QDNLS) 
is nonintegrable and cannot be analysed by the QISM. 
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2 The number state method of analysis 



In our study of these three quantum lattices we take advantage of the fact that 
each Hamihonian operator commutes with a number operator iV, which counts the 
number of fundamental quanta (bosons, g-deformed bosons, or fermions) associated 
with each degree of freedom |2^, |^] . Thus a general eigenfunction of N can be written 
in the form 

|^„)=^Q|</.i) , (2.1) 

1=1 

where p = p{n, /) is the number of different ways that n quanta can be arranged on 
/ freedoms, 10;) is the number state corresponding to a particular arrangement, and 
the {ci} are a set of p arbitrary complex constants. For example if we are placing 
two bosons (n = 2) on two freedoms {f — 2), there are three possible number states 
(P = 3): 

= |2)|0), 102) = |1)|1), and l^a) = |0)|2), which for typographical conve- 
nience we write as [20], [11], and [02]. Requiring {ipn) to satisfy the time independent 
Schrodinger equation 

H\4>„) = El^Pr.) (2.2) 
generates a px p matrix equation for the {q}'s, which determin es t he energy eigen- 



values, and the eigenvectors — upon substitution into Equation {2A) — yield the cor- 
responding stationary state wave functions. We refer to this technique as the "num- 
ber state method" (NSM) to differentiate it from the quantum inverse scattering 
method (QISM) over which it has some theoretical and computational advantages 

§. 

In effect the NSM diagonalizes the infinite Hamiltonian matrix into finite blocks 
of size p. For the bosons of Hi and the q-deformed bosons of H2 

{n + f-iy. 

P'=P'= n!(/-l)! ■ (2.3) 
Only one of the fermions of Ht, can be placed on a single freedom so in this case 

n!(/ - n)\ 

A moment's reflection will convince the reader that p can easily become inconve- 
niently large, but this is not a problem particular to the number state method. 
Exact quantum wave functions on lattices are complicated objects, and the same 
difficulty appears with the QISM [^. Each of our three systems has translational 
symmetry so these p x p blocks can be further diagonalized into smaller blocks 
with fixed values of crystal momentum fc, where r = exp (ifc) is an eigenvalue of 
the translation operator T. T is defined by the property T&j = ^j+i'T so that 

T[nin2 ••■"/] = [nfrii . . . 71/-1]. 

In this manner we find E E{k) for every allowed momentum state with a 
minimum of computational effort. 

As a simple example of this method, consider the first quantum level (n — 1). 
The energy bands for all three models are identical because there are just / ways 
that a single quantum can be placed on / freedoms. Thus we find 

|^i(ft)) = ^^(e"=r)'"'[100...0], (2.5) 
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and 

Ei{k) = ~2coak , (2.6) 

where k = 2niy/f and v = 0, ±1, ±2, . . . , ±(//2-l), //2 for / even and 0, ±1, ±2, . . . , ±(/- 
l)/2 for / odd. 

Since the effective mass (m*) is defined as 



we have 



E,{k)^E,{0) + ^ + O{k^), 



m* = 1/2. 



(2.7) 
(2.8) 



3 The second quantum level 

Here we show how to compute the exact eigenstates and energy eigenvalues of the 
sohton bands at the second quantum level (n = 2). To this end we display number 
operators, iVj , i = 1, 2, 3 that commute with each other and the Hamiltonian, and 
we construct the most general eigenf unctions of Ni and the translation operator T 
as a sum of products of elementary number states [|[ . 

(i) QDNLS 



Hi, the Hamiltonian of Equation (1.1), commutes with the number operator 



(3.1) 



which has eigenvalue n. As discussed in the previous section, the most general 
eigenfunction of iVi is a sum of products of elementary number states of the form 
[nin2 ■ ■ - Uf], where n = ni + n2 + ■ ■ ■ + nj. 

For n — 2 and / odd a general eigenfunction of both iVi and T is 



IV'2) = ^{ciE,ii(rr)^-M20---o] + c2E,ii(rr)^-Miio---o]+ 

+C3 EJ=i (rf )J"1 [1010 • • • 0] + • • • + (3-2) 

+C(/+i)/2 EW^rfy^'i^O • • • 010 • • • 00]} , 

where t = exp(ifc) is the eigenvalue of T that corresponds to the wavenumber 
k = 27rij///. To ensure that (V'2|'02) = 1, it is necessary that the q's be normalized 
as 

(/+l)/2 



E 

1=1 



1 . 



(3.3) 



Requiring that -ffilV'2) = E\ip2) leads to the matrix equation Qi{t)c = Ec, where 
c = col(ci, C2, . . . , C(y+i)/2) and Q{t) is the [(/ + l)/2] x [(/ + l)/2] matrix 



/ 7 q*V2 
qV2 q* 



L(r) 



q q* 
q p 



(3.4) 
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and 

q = l + T, p=(r(-^+l)/2+r(/-l)/2). (3.5) 

For / even (which is a special case of the condition / mod n — 0), the wave 
function depends on whether the integer v in the translational eigenvalue r = 
exp(27rij///) is even or odd. For v even 



\^2) = 47{ciE,ii(^rrM20---o] + c2Ej=i(^rrMiio---o] 



Vf 



+c3Ei=iW"'[ioio 



(3.6) 



i)/2/%/2)Ei=i(rT)^-MlO- 



•010- ••00]}, 



whereas for v odd the last term in the sum is omitted. 

For / and v both even, Qi(t) is the (//2 + 1) x (//2 + 1) matrix 

q*V2 



Qi{r) 



( V 



V 








g\/2 



/ 



(3.7) 



where q is defined in Equation (3.5). For odd Qi{t) is the (//2) x (//2) matrix 
obtained from Equation (3.7) by omitting the last row and column. 

Finding the eigenvalues of tridiagonal matrices of this sort is an interesting 
exercise in analysis: full details are given Exact results can be obtained in the 
limit / — !■ oo; for finite but large / analtic correction factors can be found. 

In the limit / — )■ oo the soliton band has the energy 



E^ik) = -V7' + 16cos2 (fc/2) 



(3.8) 



Defining the binding energy, i?b, as the difference between the energy of the soliton 
band at fc = and the bottom of the continuum band, one finds that 



Eh = Vl^ + 16 - 4 , 
and with effective mass defined as in Equation ( p/zj ) 



m 



and maximum group velocity 



dE 
'dk 



k = TT/2 



(3.9) 



(3.10) 



(3.11) 



In Figure l(i) we display the energy eigenvalues calculated from Equation (3.4) 
with 7 = 3 in the limit / — > oo. The existence of a quasi-continuum and a lower 
discrete band is clear. 



(ii) QAL 



The Hamiltonian H2 commutes with the translation operator and with the num- 
ber operator ||l|, H 



No 



' In (1 + Iblb, 



(3.12) 
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(i) QDST 



(ii) QA-L 




(iii) Fennionic model 

Figure 1: Soliton bands in the three models considered in the text 



Taking the eigenvalue n of N2 to be 2 the wave functions are as in Equations (3.2) 
and (3.6). For / odd the energy eigenvalues are found from the [(/+l)/2] x [(/+l)/2] 
matrix 



Q{r) 



( 9*72 + 7/2 

9^2 + 7/2 q* 

q 



V 



q 

q p J 



(3.13) 



where q^p are as defined in ( ^.5| ). 

For / and v both even, Q2{t) is the (//2 + 1) x (//2 + 1) matrix 



Q{r) = - 







gV2 + 7/2 
9^2 + 7/2 

9 9* 



V 



9 q 

q 0/ 



(3.14) 



whereas for v odd it is the (//2) x (//2) matrix obtained from Equation ( 3.14 ) by 
omitting the final row and column. 

In the QAL case we find two soliton bands, which are shown in Figure l(ii). The 
top band corresponds to the classical AL soliton that alternates in sign between each 
lattice point. 

In the limit / ^ 00 the soliton band has energy 



2cos(fc/2)(7 + 4) 



V27 + 4 

Defining the binding energy, i?b, as above, one finds that 

2(7 + 4) 
i^b = U 1 - 4 



V27 + 4 



and the effective masses are 



m 



± 



2V2YT4 

7 + 4 



(3.15) 



(3.16) 



(3.17) 



where the "+" corresponds to the lower band and the "— " sign to the upper band. 
In this case the maximum group velocity is 



dE 



dk 



7 + 4 
\/^7 + 4 



(3.18) 



(iii) FP 



The Hamiltonian H3 commutes with the number operator 

/ 



(3.19) 
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Because of the anti-commutation relations, the order in which one inserts the 
fermions into the chain is important |Q . We define the normal ordering to be that 
in which the fermions are inserted from left to right. Hence a2a|[00] = [11], but 



4 4 [00] 



[11]. The translation operator is defined by Ta 



a^+iT, so, for 

example, T[1001] = —[1100]. As was pointed out by Dirac, this makes hand calcu- 
lations tedious, but fortunately an algebraic manipulation system like Mathematica 
can be programmed to take such sign changes into account. 
Again assuming the eigenvalue of to be 2 and / odd, a general eigenfunction 
of N3 and f is 



-^{ci (^rp-' [110 • • • 0] -I- C2 {rfy-'mo • • ■ 0] 



/ 



(3.20) 



-i)/2Ej=iW"Mio---oio---oo]} 



Note that this wavefunction has one less term than that given in Equation (3.2) 
because no more than one fermion can be assigned to any one freedom. Requiring 
Hslips) = E\ip3) leads to the matrix equation Q3(t)c = Ec, where Q3{t) is the 
[(/ - l)/2] X [(/ - l)/2] matrix 



/7 

q 



V 



(3.21) 



and q,p are as defined in (3.5). 

117 > 2, the resulting eigenvalue plot appears as in Figure l(i), but for < 7 < 2, 
the central part of the soliton band merges with the quasi-continuum as is shown 
in Figure (iii). 

For / and v both even, 



1^2) - ^{ci EU^^Ty^'^no ■■■o] + c2 E,ii(^r)^-i[ioio • • • o] 



(3.22) 



+ --- + (c(/-i)/2/V2)E,Ci(tT) 
and Q3(t) is the (//2 - 1) x (//2 - 1) matrix 



[10- • -010 •••00]} 



Q3(t) = 



/7 q* 

q 

q 



V 



q^ 




(3.23) 



0/ 



For V odd \'ijj 2) is as in Equation ( 3.22 ) with the last term omitted, and Qsir) is as 
in Equation (3.23) but with the last row and column omitted. 
In the limit f 00 the soliton band has energy M 



E{k) 





4 2 




7^ 


— cos 




7 





for 7>2cos(fc/2). 



For 7 > 2 the binding energy and effective mass are 

4 

7 



= 7- 



(3.24) 



(3.25) 
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2 ' 



(3.26) 



and the maximum group velocity is 



'dE' 
'dk 



(3.27) 



4 Asymptotic expressions for arbitrary quantum 
levels 

Proceeding as in the previous section it is possible — in principle — to construct block 
diagonalized Hamiltonian matrices for any value of the quantum number n. As was 
noted in Section 2, however, the sizes of these blocks may grow beyond the limits 
of computational convenience (or even possibility) so it is of interest to consider 
approximate calculations that are useful in asymptotic limits. In this section we 
present results for arbitrary values of the quantum level that are asymptotically 
correct for small or large values of the anharmonicity parameter. 



4.1 Small 7 

Binding energy has been studied in detail for the QDNLS and QAL systems in 
reference and with 7 <C 1 these results are identical for QDNLS and QAL. To 
calculate the binding energy in this limit it is necessary to know whether the size 
of the classical soliton is large or small compared with the number of freedoms. If 
1 ^ 7 > 24/((n + 1)/), the size of the classical soliton is smaller than / and 
both the QDNLS and the QAL equations are well approximated by the continuum 
nonlinear Schrodinger equation for which the binding energy of a quantum soliton 
is 

E^^^nin'-l). (4.1) 

To calculate the effective masses for QDNLS and QAL in the limit of small 7 
we refer to Equation ( ^.6| ) , which gives the energy as a function of the wave number 
for a single quantum. Since the classical problem is linear in this limit, the energy 
for n quanta with wave numbers: ki, k2, ■ ■ ■ , kn is just the sum: —2 cos kj. For 
the n-quantum wave function k — kj, and the lowest value of energy is found 
for fci = k2 = ■ ■ ■ ~ kn', thus in this limit 

En{k) = -2ncos(fc/7i) (4.2) 

so 

71 

m* = -, (4.3) 
and the maximum group velocity of a wave packet is 

'dE 



dk 



k=ir/2 



= 2sm(|-). (4.4) 



For the FP system the soliton band merges into the continuum at small 7 so i?b 
and m* are not defined in this limit. 



4.2 Large 7 

In the limit of large 7 it is evident that for n = 2 the dominant elements of the 
matrices Qi, Q2, and Q3, which are displayed in Equations (3.4), (3.12), and (3.21), 
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are the 2x2 submatrices in the upper left hand corners. For arbitrary values of n, 
perturbation theory in 7"^ shows that to calculate the leading fc-dependent terms 
it is only necessary to consider the sequence of interactions 

[n] ^ [n - 1, 1] ^ [n - 2, 2] 4-> • • • <-> [2, n - 2] 4^ [1, n ~ I] ^ [n] . (4.5) 

For the QDNLS and QAL systems this sequence requires an approximate wave 
function of the form 

\^n) = ^{ciEj=iW-M'^o-..o] + c2E,ii(^^p-M("-i)io---o]+ 
+C3 [(" - 1)00 . . . 1] + • • • + 

+c„E,ii(^rrMK2)W2)o. • -0]} , 

(4.6) 

with n even, and 

\^n) = ^{ciEj=iW-M'^o...o] + c2E,ii(^r)^'-M("-i)io---o]+ 
+C3 [(" - 1)00 . . . 1] + • • • + 

+c„-i Ej=i(^T')^'"M((^i + l)/2)((n - l)/2)0 • • - 0] , + 
+c„E,ii(^r)^"M((" + l)/2)0. • •0((n- l)/2)]} 

(4.7) 

for n odd. 

(i) QDNLS 

As an illustrative example, consider the QDNLS equation with n even for which 
the n X n matrix — Qi(t) formed from the approximate wave function in Equation 



( |4.6| ) is 



n(n— 1) 
2 










„,(n-l)(«-2) 
' 2 












„(«-l)(«-2) 
' 2 


V3(n-2) 



7( x-t + l) T-Vt(t + 1) 

V vW+1) ^yt(f+l) 7(^-§) / 

(4.8) 

We have no exact results for the eigenvalues of such approximate matrices for 
arbitrary values of n. However we have done extensive investigations of the asymp- 
totic expansions of the eigenvalues in powers of I/7 for finite values of n using 
Mathematica |^ and perturbation theory. Clearly the most negative eigenvalue 
has a leading term —n(n — l)/27. Higher order terms can be calculated but they 
are of secondary interest except for the lowest order term that depends on k. Per- 
turbation theory shows that this term is of order 1/7"^^. In summary, the most 
negative root of the determinental equation, det \Q\ — 1E\ = 0, defines a soliton 
band of the form 

En{k) = --n{n - 1)7 - _ j cosfc , (4.9) 
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where the symbol indicates that the first term on the right hand side is correct 
to 0(7"^) while the second (fc-dependent) term is correct to 0(7-'^""'). In other 
words, terms of order 0(1/7) that do not depend on k have been dropped. The 
second term on the right comes from the perturbation correction. 
Equation (|4.9|) implies that for the QDNLS equation 



Eh = ^n{n - 1)7 , 
2n 



and 



'dE' 

'dk 



2n 



-./2 ("-1)!7"-^' 



(4.10) 
(4.11) 

(4.12) 



(ii) QAL 

For the QAL equation the wave functions are of the forms given in Equations 
( |4.6| ) and (4.7) except that it is convenient to permute the order of terms such 
that {ci, C2, 



J Cn } ^ { _ 3 , C< 



n — 5 1 



. ,C3,ci,C2,C4,...,c„_2,c„,c„_i}. With this 



rearrangement the corresponding n x n matrix is to leading order in 7 



Q2{r) 





p 












P 


\ 







p 















p 





p 





















p 

















p 





PT- 


1 














Pr 





/ 



(4.13) 



where 



Pin) . (I) 



(n-l)/2 



(4.14) 



Now define Qn ~ det[/3 ^Q2{t) — el]- Expanding Qn by final rows and columns 
gives eventually that 

Qn = (-l)"+i2cosfc-2P„_2 + eP„_i 

where Pn-i is the tridiagonal determinant formed by removing the final row and 
column of Q„. 

A standard calculation shows that P,i(e) is a polynomial satisfying the recursion 
relation 

Pn = ePn-l - Pn^2 , (4.15) 

where Pi = e and P2 = — Thus P„(e) = 2C/„(e/2), where Un{x) is a Chebyshev 
polynomial of the second kind. Hence 



Qn = (-l)"+'2cOsfc + C/„-C/„_2 

= (-l)"+i2cosfc + 2T„(|) 
Where T„(a;) is a Chebyshev polynomial of the first kind. 



(4.16) 



A further short calculation shows that the eigenvalues of (4.13) are 



Enik) 



-2cos 



7\ ("-i)/2 



(i) 



(4.17) 
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so for n > 1 



and 



7N ("-i)/2 



„2 /2X("-l)/2 



m = 



(4.18) 



(4.19) 



2 V7. 

The maximum group velocity on the QAL soUton band occurs at the band edge 
{k = it), therefore 



dE 
Ik 



7r\ /7\("-i)/2 



fc=7r 



(4.20) 



(ill) FP 

For the FP model a translationally invariant wave function in the large 7 limit 



IS 



l^n) = ^{ciEtiW'-Mii---io---o]+ 

+C2 ELi (tT )J-Mll • • • 1010 • • • 0] + 



(4.21) 



+C3 Ei=i (tT)^-! [Oil • • • 10 • • • 01]+ 
+ • • • + c„ EJ=i W"'[li;;_l li;;_l • • • 0]} , 

7i/2 n/2 

for n even. For n odd there will be a corresponding expression with the last two 
terms given by 

cn-i EJ=i W"' [ll;;JO il^O • • • 0]+ 

(n-l)/2 (n+l)/2 

cn E J=i [rfy-^ [n^ 11_^ • • • 0] } . 

(n+l)/2 (n-l)/1 

Again it is convenient to reorder to get a periodic tridiagonal matrix. For n even, 
take {ci,C2, . . . ,c„} {ci, C3, C5 . . . , c„/2, c„/2_i, c„/2_3, • • • , C2}, with a similar 
permutation in the n odd case. This leads to a n x n matrix ~Qz{t) 



/(n-l)7 

T (n-2)7 1 
1 (n-2)7 1 







1 \ 




V 







1 (n - 2)7 1 
1 (71-2)7/ 



(4.22) 



The characteristic equation is calculated in a similar way to the QAL case. If 
we define Qn = Aet{—Q^{T) + (n — 2)7 — el], then 



Qn = ^ (7 + e)f/„_i(e/2) - 2f/„_2(e/2) - (-l)"2cosA; = 



12 



To solve this polynomial equation for e by a perturbation series around j — oo, put 

7 — e^^, z = — ee, to get 

M 

- 1) + (-l)"+ie"2cosfc + e2'"z2(M-m)p^^^) ^ q 

m— 1 

where M = [n/2j and the Pm are linear functions of z with constant coefficients. 
A series expansion shows that one root is z = 1 + O(e^) and the others are 0{e). 
If we seek t he lo west order A:— dependent correction to the z — 1 root, along the 
lines of eqn (O), it is not difficult to show that this occurs at 0(e") and is equal 
to e"2cosfc. Resubstituting this back into the original problem, we have that the 
soliton band is given by 

£;„(fc) = -(n- 1)7- -^cosfc, (4.23) 

where the symbol "=" has the same meaning as that discussed following (J^^ . This 
implies that 

Eb = {n- 1)7 , (4.24) 
m* = -^Y- ' and (4.25) 



dE 



dk 



— 1 

fe=7r/2 ' 



(4.26) 



5 Soliton wave packets 

The picture that emerges from our studies of one dimensional quantum lattices 
with / degrees of freedom and translational symmetry is as follows. For each value 
of the principle quantum number n and wave number k, there is a lowest energy 
eigenvalue. These / lowest eigenvalues lie on a band (see Figure 1) 

E = Enik), (5.1) 

where the wave number 

k='-f, (5.2) 

and 1/ = 0, ±1, ±2, . . . , ±(//2 - 1), //2 for / even and 0, ±1, ±2, . . . , ±(/ - l)/2 for 
/ odd. Each energy eigenvalue corresponds to a pure eigenstate \ipn{k)), which is 
normalized as (V'n(^)IV'n(fc)) = 1- 

For H being one of the Hamiltonian operators in Equations (p~4|), ( |1.2D or ( |l.3D , 
solutions of the time dependent Schrodinger equation 

il\^{t)) = H\^{t)) (5.3) 

can be constructed as sums over the quantum number and the wave number. Thus 
for QDNLS and QAL 

ll'(t)) =^a„^G„(fc)|^„(fc))exp(-ii;„(fc)i), (5.4) 

n k 

where n = 0, 1, 2, . . . , oo, and k takes the values between — tt to vr that are indicated 
in Equation (5.2). Since (4'(t)|5'(i)) = 1, both the a„ and the G„(fc) are sets of 
complex numbers that satisfy the normalization conditions 

^|a„|2 = l and |G„(fc)|2 = 1 . (5.5) 
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It should be noted that Equation (5.4) does not represent the most general wave 



function that satisfies Equation (5.3) because it is constructed only from eigenstates 
with eigenvalues on the soliton bands. It is seen from Figure 1 that the system has 
many other eigenstates, and these are excluded from |^'(i)) as defined in Equation 



(5.4), which is characterized by two interdependent properties: 



• For given values of n, fc, and the expansion coefficients, |^(i)) has the lowest 
energy, and 

• Under the same conditions, \'^{t)) has the highest probability of quanta being 
located near each other. 

These properties are the basis for referring to l'I'(i)) as a "soliton wave packet". 



An effect of dispersion in the wave packet of Equation (5.4) is to introduce 
uncertainties in position (j) and crystal momentum (fc) that satisfy the Heisenberg 
relation. Such uncertainties appear in a natural way because G„(fc) is a discrete 
Fourier transform of the soliton pulse shape so 

AfcxAj-1. (5.6) 

This point has been neglected at times in discussions of "Davydov's soliton" , which 
has been proposed as a polaronic means for transport of energy or charge in protein 

In this section we present some properties of soliton wave packets. In addition 
to considering the differences between the QDNLS, QAL, and FP models, we also 
distinguish between spectral problems, where n has a small, fixed value, and the cor- 
respondence limit for QDNLS and QAL, where values of n are large compared with 
unity and |5'(i)) is approximated by the solution of a classical nonlinear equation. 

5.1 Spectral problems 

The ground state (n = 0) energy is zero for all three models, and the corresponding 
eigenstate is j-^o) = [000 • • • 0]. For n — 1 the wave function of an exact eigenstate 
is given in Equation ( ^.5[), again for all three models, and the energy eigenvalues 
are given in Equation ( ^.61 ). For n = 2 exact expressions for \'>p2{k)) and E2{k) are 
presented in Section 3, and in Section 4 we have derived approximate expressions 
for l^ri(fc)) and En{k) that hold in the limits 7 <C 1 and 7^1. 



When the expression for |^(f)) in Equation (5.4) represents the wave function 
at a particular energy level (say n — rip), then 

In measurements involving infra-red or Raman spectra on molecular crystals, the 
wave length of the interacting radiation is much larger than the dimensions of a 
unit cell so it is often reasonable to assume fc « 0, but one can arrange experiments 
for which this is not the case. 

The results of Section 3 permit exact calculations of features for the three tran- 
sitions: 

l^o) l^i), 
l^o) ^ 1^2) , 

IV'l)^l^2). 

From the results of Section 4, which obtain in the limit 7 3> 1, approximate features 
can be calculated for the transitions: 
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■ ■ ■ etc. • • • 

\lpm) \tpn) ■ 

5.2 Large wave packets 

(i) QDNLS solitons in the Hartree approximation 

In the continuum limit this problem has been studied in detail by Lai and Haus 
and by Wright as a model for the propagation of solitons on an optical 
fiber. The exact solution exhibits two quantum effects: phase spreading, which 
is caused by different values of n in the wave packet, and dispersion, caused by 
different values of k. Lai and Haus have shown that the effects of dispersion become 
negligible compared with phase spreading at large values of the average number of 
quanta (bosons) in the soliton. 

For the QDNLS equation with n 3> 1, the soliton wave function at a particular 
value of n 

l^nit)) = G„(fc)|^„(fc))exp[-i£;„(fc)i] (5.8) 
fe 

is close to the Hartree approximation |2£| 



f 



|0) , (5.9) 



where $„j (i) is a solution of 



- +($„j+i + $„j_i)+7(n-l)|$„,jf $„,j =0 (5.10) 

that satisfies the normalization condition 

/ 

El*-^WI' = 1- (5.11) 
i=i 

Since ''^^'4'n{j,t)\bj\^pnij,t))^-^'> = 0, it is convenient to choose the {a„} in 
Equation (^.4|) so that |^(i)) is the coherent wave packet 



^it)) = E ^exp (-^) \Mj,t)r^ . (5.12) 



^/2 
) 

In the correspondence limit, no S> 1, 

{^'it)\b,\^it)) ^Mt), (5.13) 
where Aj (t) is a solution of the classical equation 

i^ + {A,+,+A,_,)+j\A,\'A, = (5.14) 

with normalization 

/ 

J2\A,\^=no. (5.15) 
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For 7 > 0(1) both classical solutions of Equation ( ^.14 ) and solutions of Equa- 
tion ( 5.101 ) do not propagate; they are pinned to the lattice by the Peierls barrier 
The exact expression for |4'(i)) is seen from Equation (4.12) to have a maximum 
wave packet velocity 

(5.16) 



(no - l)!7"o-i 



for riQ ^ 1. Thus the exact quantum mechanical solution also becomes pinned in 
the correspondence limit in the sense that Vi„ — > strongly as np — + cxd. 

(ii) QAL solitons 

For 7 <C 1 the QAL equation approaches the QDNLS equation and the wave 
function of a QAL soliton is close to that given by Equations (5.E) and ( 5.12| ). 
For 7 > 0(1) there is — to our knowledge — no Ha rtree appr oximation, and the two 
systems are quite different. From Equations (L€) and (4.7) 



and 



= py'Ty^' MO • • • 0] + 0(7-^) 

6jfJ-^[nOO---0] = (3{n)f^-^[{n- 1)00 ■■■0], 



(5.17) 



(5.18) 



where (3{n) is defined in Equation (4.14). Thus it is possible to construct 



where \'^j{t)) is the coherent state 



n(n-l)/4 



$"(t)T^-i[nOO---0], 



for which 
and 



b,\^,{t))^^,it)\^,it)) 



CO f 

EE 

n=0 j = l 

Thus the Heisenberg operator equation 
db 



(3{n) 



-1/2 



i^ + (6,+i+6,_i)(l + |6]6,)-0. 



and Equation ( ^.21 ) imply that 



dt 



(a>,+i + a>,_i)(i + Jci>,f) = o. 



(5.19) 



(5.20) 
(5.21) 
(5.22) 

(5.23) 
(5.24) 



From Equation ([4.20|) we see that the maximum group velo city is unbounded 
as n ^ 00. This is consistent with the fact that Equation (5.24) supports moving 
solitons that do not become pinned as 7 increases m. 



(iii) FP solitons 
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Having arisen in the context of classical dynamics, the term "soliton" is es- 
sentially a classical concept, and a "quantum soliton" is often considered to be an 
object that becomes a classical soliton in the correspondence limit. For the fermionic 
polaron model, however, there is no correspondence limit because the number of 
fermions is at most equal to the number of freedoms. Nonetheless one can construct 
a wave function of the form of Equation (5.4), where n — 0, 1, 2, . . . , /. Since this 
wave packet shares many properties with Equation (5.4), and the FP equation can 
be analyzed using the q uan tum inverse scattering method | pO| , it is appropriate to 
call in Equation (5^) a quantum soliton. 

Since there is no correspondence limit for an FP soliton, it is not in app ropriate 
to choose a particular value for n (say np) for the sum in Equation (^J). Thus 
On = <5Ti,no; and for 7 ^ 1 an approximate picture of our soliton (with velocity v) 
becomes 

no 

[0000 • • • 0000 yC~1 0000 •• • 0000] . 



However this diagram does not show uncertainties in position (j) and momentum 
(fc), which are required by the Heisenberg principle and included in the structure 
of Equation (5.4). For n ~ no this equation takes the form 



\'^no{t)) = ^G„o(fc)l^«o(fc))exp[-i£;„o(fc)i] 



(5.25) 



where, from Equation (4.21), 



|V^„„(fc)) = -^^(c'^Tri[ll---100---0] 



and, from Equation ( 4.23 ), 



-(no - 1)7 - 



v'lo - 1 



cosfc . 



(5.26) 



(5.27) 



Since G„(,(fc) is the probability amplitude for the FP soliton to have momentum /c, 
Equation ( ^.25| ) can be written as 



where 



{t)) = l=Y,Fn„ (j, t)f ^-'[11 • • • 100 • • • 0] , 

Fnoij, t) = J2 Gno(fc)exp[i(fc(j - 1) - E„,{k)t)] . 



(5.28) 



(5.29) 



Thus uncertainties in position and momentum are related by Equation (5.6). If k 
lies within the range 

Ak 

fc = fco ± — , (5.30) 



then the soliton speed is in the range 
Av 2 



v± 



1 



"0-1 



Ak , . Ak 

sinfc cos ± cosK sin 

2 2 



(5.31) 



Since wave functions of an FP soliton are dominated by components of the form 

[•••OOIT^^^OO---], 

n times 
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it becomes more extended as the number of fermions increases. For an infinite num- 
ber of freedoms on the lattice, one could allow the number of fermions in Equation 



(5.4) to grow without bound, but the size of the soliton would also become infinite. 

For n = f, \ipf) — [111 .. .1] with energy Ef = —7/, while for f — 0, \4>q) — 
[000 ... 0] with energy Eq = 0. In general the exact eigenfunctions for n = f — m will 
have the same structure as those for n = m, but with the I's and O's interchanged. 
Thus for n = f — m, where 771 = 0, 1, ... , /, the corresponding energy will be 

Ef^mik) = -lU - 2777) + E„,{k) , (5.32) 

where Em{k) is the soliton energy band in the case n = m. 



6 Conclusions 

The primary aim of this work is to make clear the concept of a quantum soliton 
by presenting several specific examples. Lattice solitons are of particular interest in 
this effort because it is possible — in certain cases — to examine details of the corre- 
sponding wave functions and to appreciate their complex character. The properties 
of quantum solitons are found to depend strongly upon the level of anharmonicity 
and the commutation relations that characterize the fundamental quanta. 

The number state method of analysis, which has been somewhat obscured by 
the quantum inverse scattering method, is found to be a useful computational tool 
for such problems as well as a helpful theoretical perspective. 

Whereas previous publications have concentrated on solving the quantum prob- 
lem with only a few quanta using the number state method, we have shown that 
an approximate version of this method, valid in the limit of large nonlinearity, is 
also useful for larger values of 77. In particular, we have found general expressions 
for the binding energy, the effective mass and the maximum group velocity in this 
limit for arbitrary 77 > 1. 



model 


Eb 


777* 


Kn 


QDNLS 


^n{n — 1)7 


2n 1 


2n 

(n-l)!7"-i 


QAL 


2(^)("-i)/2 


„ / X (n-l)/2 
2 {1) 


n ^71' \ 2 / 


FP 


(77 - 1)7 


2 


2 

7"-^ 



We note that the results for the group velocities confirm that the QDNLS and 
the FP solitons get pinned for large 7 whereas this is not the case for QAL solitons, 
in agreement with the findings of classical soliton theory. 

This method can also be applied to Hubbard models jlSj, which are of inter- 
est as theories of superconductivity and ferromagnetism. This application will be 
discussed elsewhere. 

Finally we note that the band structures shown in Figure 1 may be somewhat 
misleading because they are for two fundamental quanta (77 = 2). With 77 > 3, 
exploratory numerical and theoretical studies show that additional bands reside in 
the gap between the lowest (soliton) band and the principal continuum band. The 
significance of these bands will be the subject of future research. 
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